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Abstract 

We show that the pattern of tunneUing rates can display a vivid and regular pattern when the classical dynamics 
is of mixed chaotic/regular type. We consider the situation in which the dominant tunnelling route connects to a 
stable periodic orbit and this orbit is surrounded by a regular island which supports a number of quantum states. 
We derive an explicit semiclassical expression for the positions and tunnelling rates of these states by use of a 
complexified trace formula. 
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Tunnelling in systems whose classical limit dis- 
plays a mixture of chaotic and integrable be- 
haviour ||l|-Q] is often quite complex and impos- 
sible to predict analytically. Much attention has 
been paid recently, for example, to the regime of 
chaos- assisted tunnelling ||] in which dynamical 
tunnelling occurs between quasimodes supported 
in integrable island-chains embedded in a chaotic 
sea. By contrast, we report on a remarkably or- 
dered structure that appears in the tunnelling 
behaviour of a particular kind of mixed system 
and give analytical estimates for the correspond- 
ing tunnelling rates. It is distinct from the case 
of chaos-assisted tunnelling because tunnelling is 
through an energetic barrier rather than through 
dynamical barriers such as KAM tori. The spe- 
cial feature of these systems is that the complex 
orbit which defines the optimal tunnelling route 
across the barrier connects to a stable periodic 
orbit at the centre of an island chain (which is 
generally embedded in a chaotic sea). 

The ordered nature of these tunnelling rates is 
immediately evident in Fig. m where we show the 
numerically obtained splittings between quasi- 
doublets of the double- well potential V{x,y) ~ 
(x^ - 1)4 + x2y2 + 2y2/5. We have held the en- 
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Figure 1. Each dot represents a splitting doublet. 
The horizontal coordinate is the mean value of g = 
1/h and the vertical coordinate is the splitting. 

ergy fixed (at E = 9/10) and found the values 
of g = 1/h for which this is an energy level. 
The resulting spectrum of q values is equivalent 
in most respects to a standard energy-spectrum, 
with the advantage that the classical dynamics is 
fixed throughout. 

The largest splittings in Fig. |l| are highly or- 
dered, forming a regular progression of families 
which grow larger in number higher in the spec- 
trum. These correspond to states supported near 
the centre of the island and we will offer a sim- 
ple analytical prediction for them. The small- 
est splittings in Fig. n^ form a disordered jum- 



ble. These correspond to states supported in the 
chaotic sea and dynamicaUy excluded from the 
main tunneUing route. 

To analyse the ordered sequence we use a 
method developed in |^,g| and used until now pri- 
marily to understand predominantly chaotic sys- 
tems. We first present the analysis for the case 
that q = 1/h is held fixed and an energy-spectrum 
is computed. We then give the simple extension 
for the fixed energy g-spectrum, appropriate for 
the results of Fig. ffl. We call the mean energy of 
the n'th doublet E^ and the corresponding split- 
ting AEn and define the following dimensionless 
spectral function 



/(-B,g)=^A£;„5(^-K). 
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There is an analogous definition for metastable 
wells which have extremely narrow resonances, in 
which the widths r„ play the role of splittings. 
Such a system was studied in for a fully in- 
tegrable system; the structure of the spectrum 
was like that shown here except with no irregular 
jumble at the bottom. While similar in outline, 
the detailed method of analysis was rather dif- 
ferent, making use of the action-angle variables 
which exist in that situation. 

In pM| we approximate M) semiclassically as 
a sum over complex tunnelling orbits which tra- 
verse the barrier (in analogy to Gutzwiller's for- 
mula for the density of states using real orbits [|j). 
We shall consider the special case in which there 
is an additional reflection symmetry such that the 
dominant tunnelling route lies on the symmetry 
axis so that it connects smoothly to a real periodic 
orbit. We then identify three distinct contribu- 
tions to f{E,q). There is the so-called instanton 
which has a purely imaginary action iKq, lives 
under the barrier and runs along the symmetry 
axis between the classical turning points. It is 
important for determining the mean behaviour of 
the splittings but does not affect the fluctuation 
effects which we are trying to capture here. The 
second contribution comes from orbits which ex- 
ecute real dynamics along the real periodic orbit 
lying on the symmetry axis in addition to the 
instanton dynamics beneath the well. We imag- 
ine an orbit which starts at one of the turning 



points, executes r repetitions of the real periodic 
orbit and then tunnels along the instanton path 
to finish at the other turning point. This orbit 
has a complex action S = rSo + iKq where Sq 
is the real action of the real periodic orbit. The 
contribution to f{E,q) is given by H 
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The matrices Wq and Mq are the monodromy 
matrices of the instanton and of the real peri- 
odic orbit respectively; the composite orbit has 
a monodromy matrix which is simply a product 
of these. The third contribution, discussed in |^, 
comes from homoclinic orbits which explore the 
real wells far away from the symmetry axis. They 
play no role in the present discussion. 

For fully developed chaos, all periodic orbits 
are unstable. The denominator of (0) then decays 
exponentially with r and large repetitions are nu- 
merically unimportant. If the orbit is stable, how- 
ever, there is no corresponding suppression. Sin- 
gularities corresponding to distinct states arise 
when the expression is summed. This follows 
very closely the analogous development of Miller 
and Voros |9|-|ll[ for the Gutzwiller trace formula 
when there is a stable orbit, except here we find 
splittings in addition to the positions of energies. 
(As was done in one dimension by Miller [[l2[.) 

In the stable case, Mq has eigenvalues e**" 
on the unit circle. (In higher dimensions, there 
would be a number of such eigenvalues and the 
theory would be generalised accordingly.) Let the 
diagonal matrix elements of Wq in the complex 
eigenbasis of Mq be A and B so that 



det{WQM[^ - I) 



= TrWoM^-2 

= Ae*™ + Be"'™ - 2. (3) 



Since the instanton's period is imaginary, com- 
plex conjugation acts as a time-reversal opera- 
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We therefore 



tion and we find that Wq 
conclude that (TrVFoM^)* = Tt{WoMq'')-^ = 
TiWqMq'' (the latter equahty holds because ev- 
ery symplectic matrix is conjugate to its inverse). 
Comparing this with (||) we conclude that A and 
B are real; we discuss how they are computed and 
offer a geometrical interpretation in the appendix. 



We now make use of the generating function of 
the Legendre polynomials to conclude 
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where we assume without loss of generality that 
B is the larger in magnitude of [A, B) and we let 
Qk{z) denote the polynomial 



Qk{z] 
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Using this in (0) and summing the resulting geo- 
metric series in r for each k we get, 
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$fe = qSo + {k + l/2)a. 
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Semiclassical energy levels are found when the 
distribution above has poles and are implicit so- 
lutions Emk of <i>fc = 2Trm. From the residues we 
recover estimates of the corresponding splittings. 
This is a form of torus quantisation in which k 
is a transverse quantum number, treated in har- 
monic approximation, and m counts nodes along 
the orbit. The corresponding states are localised 
on the tori surrounding the stable periodic orbit 
and we find that their respective tunnelling rates 
are much larger than those of other states. 
Note that near Emk we can write 
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where we have used that the period is Tq = 
dSo/dE. Using the standard identity 
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E — E,nk 

we conclude 
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where the notation now stresses that aj, depends 
on the transverse dynamics of the instanton and 



its real extension through the monodromy matri- 
ces (Wo,Mo). All classical quantities are evalu- 
ated at energy E^k- This formula has the same 
form as in one-dimensional tunnelling [|l2| except 
for the additional factor ak{Wo, Mq). We will find 
this factor to be of order unity when k is small 
and to decrease as k increases. 

In m we reported some specific numerical re- 
sults for the energy quantisation. We now extend 
this result to the the g-spectrum, the advantage 
being that the classical quantities Kq, Wq and Mq 
are constant. The function f{E,q) can equally 
well be interpreted as a function of q at fixed E 
so that (||) still applies. Solving for the poles and 
residues as above, we conclude 
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The sequences of Fig. |l| can now be interpreted 
in terms of the quantum numbers m and k. 

The central states with k = correspond to the 
largest splittings, about 35 times larger than the 
local average. In Fig. |l| they are the uppermost 
curve of points (along which m varies) . Keeping 
m fixed and letting k increase one gets a sequence, 
which appears as a left-sloping shoulder in Fig. ^, 
along which q and Aq decrease. 
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Figure 2. The open circles are numerically-computed 
splittings. The filled circles are the semiclassical pre- 
dictions for the family with m = 40 and < fc < 6. 

In Fig.|| we show a subset of the spectrum with 
the semiclassical predictions for m — AQ and a se- 
quence of k values (using A — —0.861, B = 4.043 
and a = 2.783.) Clearly the small k states are 



well reproduced. Our analysis essentially extrap- 
olates the properties of the central periodic orbit 
to the entire island Q, and this is less accurate 
for the large k states which are localised further 
from the periodic orbit. Reproducing the large k 
values would require a more sophisticated anal- 
ysis [|ll[, although our formalism does at least 
capture the correct qualitative behaviour of these 
states. The irregular jumble of splittings at the 
bottom of the figure corresponds to states in the 
chaotic part of phase space. No simple theory ex- 
ists for them though one could well imagine that 
the formalism of chaos-assisted tunnelling, in par- 
ticular the interplay between regular and chaotic 
states, might be of use in describing them. 

As remarked, the number of well-defined k 
states increases as we go up in the spectrum. This 
is because the number of states which the regular 
island can support increases as h decreases. Also, 
there are occasional irregularities in the lattice of 
regular states, for example the fc = state near 
q = 18. These are due to near degeneracies be- 
tween the regular state and some other state — 
either another regular one or a chaotic one. The 
actual eigenstates are then strongly mixed and 
hence so are the tunnelling rates. 

Appendix 

To apply Eq. ( pT| ) in practice one needs the 
parameters A and B; we give here simple basis- 
independent expressions for them. In particular, 
we note that it is not explicitly necessary to trans- 
form Wo to the eigenbasis of Mq. We calculate A 
and B as the smaller and larger respectively of 



A or B = cosh /3 ± 7 sinh /3 
where TrWo = 2 cosh /? and, 
ImTrW^oMo 
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This is obtained by expressing Wq = e""^^'^^ and 
Mq = e"''^ where J is the unit symplectic matrix 
and H and K are real, positive-definite, 2x2 
symmetric matrices normalised so that dct H = 
det K = 1. Expanding 



and similarly for AIq^ — e^"'-'^ , one recovers (g) 
with A and B as given above. 

The factor 7 has the following geometric inter- 
pretation. The action K -^ MKM"^ of 2 x 2 
symplectic matrices M on symmetric matrices K 
can be identified with (2-t- l)-dimensional Lorentz 
transformations (since the relevant Lie algebras 
are isomorphic), the invariant det K playing the 
role of proper time. The matrices H and K de- 
fine unit time-like (2 -I- l)-vectors X = {x,y,t) 
and S — (^,?7,t) respectively. For example 
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x^-y^ = l, (15) 



Wo = e^'^'^" = cosh/3 - iJiJ sinh /3 



and similarly for K . One then observes that 7 = 
-TtJHJK/2 = {X, E) =tT -x^- yq. This can 
be interpreted as the dilation factor to boost the 
rest-frame of X to that of S. In particular, 7 > 1. 
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